function ds = ZermeloSys(t, s, tu, u)
u = interp1(tu, u, t);
ds = [cos(u)+s(2);sin(u)];
